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Abstract: A new method to reconstruct the 3-dimensional structure of extensive air show- 
ers, seen by fluorescence detectors, is proposed. The observation of the shower is done in 
2-dimensional pixels, for consecutive time bins. Time corresponds to a third dimension. As- 
suming that the cosmic ray shower propagates as a plane wave front moving at the speed of 
light, a complex 3D volume in space can be associated to each measured charge (per pixel 
and time bin). The 3D description in space allows a simultaneous access to the longitudinal 
and lateral profiles of each shower. In the case that several eyes observe the same shower, 
the method gives a straight-forward combination of all the information. This method is in an 
early phase of development and is not used for the general reconstruction of the Auger data. 
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Introduction 

The Pierre Auger Observatory will provide a 
large set of cosmic ray data to be analysed 
in multiple perspectives, ranging from particle 
physics to cosmology. A detailed understand- 
ing of the data is crucial. 

The method proposed in this contribution aims 
at reconstructing fine-details of the individ- 
ual shower structure and being sensitive to 
both lateral and longitudinal profiles simulta- 
neously, from the fluorescence light emission. 
It is not the standard reconstruction method 
used in Auger, but is built on top of the ex- 
isting methods that provide the geometry of 
the shower line and of the longitudinal profiles 
with great accuracy. 

The Auger Fluorescence Detector (FD) is com- 
posed of four eyes, each with six telescopes, 
observing the atmosphere over the centre of 
the Surface Detector array from different direc- 
tions, with each eye covering elevations 2° to 
32° and 180° in azimuth. The data comprises 
cosmic ray events up to the highest energies, at 
very different distances from the detector eyes 
- ranging from 1 to above 30 km. 



The Auger FD data are collected in pixels of 
approximate angular dimensions of 1.5° and in 
100 ns time bins. The standard reconstruction 
is based first on the pixel angular information, 
to find the plane that contains the shower axis 
and the observing eye (the Shower-Detector 
Plane, SDP). The centroid time found for each 
pixel is then used to find the axis line within 
this plane, the minimum approach distance to 
the detector, and the reference time at which it 
occurred (Tq). The shower geometry is simply 
given by this line (details of the calculation are 
given in [1], the reconstruction can also use sev- 
eral eyes and the surface detector, as explained 
in [4]). In the above procedure pixels which do 
not observe the SDP cannot be used for the ge- 
ometry reconstruction and the time structure 
of the signal inside each pixel is neglected. Al- 
though for distant showers the line approach is 
clearly sufficient, for close-by showers relevant 
information can be lost. 

To measure the energy deposited by the pri- 
mary cosmic ray, a Gaisser-Hillas function is 
then fitted on this line. To include the effect of 
non-perfect optics and the fact that the shower 
in not only a line, the contribution of near-by 
pixels is considered at this stage, by summing 
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signals for each time bin. Some of the infor- 
mation important for the longitudinal profile 
and energy reconstruction is restored in this 
procedure. 

However, the information that can be present 
in the lateral profile of the shower, in possi- 
ble asymmetries or local fine structures is lost 
in the standard reconstruction - this is the 
main motivation to find another reconstruction 
method that takes all the available information 
into account. The method described below is 
under test in simulation (and laser data), and 
will soon be tested in real data. 

The 3D method 

Geometry Reconstruction 

The basic idea of the method is to propagate 
the available 3D information - binned in two 
angular dimensions and one time dimension - 
into 3D volumes in space. These volumes con- 
tain the points from which the observed fluo- 
rescence photons were emitted. 

The cosmic ray shower is assumed to propagate 
as a plane wave front, moving at the speed of 
light. For each direction observed in the de- 
tector, there is a line in space that in general 
does not intercept the shower axis. The ob- 
servation time, on the other hand, is a sum of 
the emission time of the photon and its prop- 
agation time to the detector. The emission 
time is characteristic of all the particles in a 
disk moving coherently along the shower axis 
at constant velocity c. 

Starting from a given hypothesis for the shower 
axis and core distance, each observation direc- 
tion at each time corresponds to a unique point 
in space. The axis and core coming from the 
standard reconstruction (being it mono, stereo 
or hybrid when available) are used as a first 
hypothesis, and, with the geometry fixed, the 
time of closest approach between the shower 
line and the detector, Tq, is found using the 
pixel centroids and directions within the SDP, 
as before. 

The position of the centre of each volume (cor- 
responding to 6j , 4>j of the pixel centre and U 




Figure 1: 3D view of a simulated event of 
2^018.5 gY gggjj 2 km from the telescope. The 
shapes of the 3D volumes are shown and the 
color code corresponds to the detection time. 
The visualisation is done with the map3d pack- 
age [3]. 



for each time bin) is then found for 



being aj the angle between the observation di- 
rection and the shower axis. 

The borders of each volume are given by the 
same procedures using the 9, (j) of each vertex 
and ti±50 ns. Notice that the Auger pixels are 
regular hexagons layed on a spherical surface 
- irregular in 9,(j) - and that the observation 
times correspond to constant emission times 
along the bisectrix between the shower axis and 
the observation line. 

As the time differences between pixels give in- 
formation on the shower axis location, the time 
duration of the signal within each pixel can 



30th International Cosmic Ray Conference 



also provide the same information. In addi- 
tion, larger time signals can also correspond 
to larger lateral dimensions of the shower. To 
profit from the extra information, the charge 
observed in each volume is used to refine the 
shower axis (and core) determination. 

This is done by calculating the "centre-of- 
mass" and "inertia axis" of the 3D object 
composed by all the volumes, considering the 
charge to be, in first approximation, propor- 
tional the number of particles, and to the 
"mass". 

The main "inertia axis" gives the shower di- 
rection, replacing the shower axis (the shower 
core being fixed by the "centre-of- mass"). The 
two other axes give a first measurement of the 
average lateral distribution, the fact that they 
are usually equal shows that there is no signif- 
icant bias in the use of time as a third dimen- 
sion. 

To a new axis and core hypothesis correspond 
new SDP and Tq (the latest to be recalculated 
from data), and so new positions and shapes 
for each volume. The procedure can thus be 
iterated until reasonable stability is obtained 
for the axis direction and core location. The 
resolutions on core location and axis directions 
obtained after the iterations are comparable 
with the ones from the standard reconstruction 
used as inputs, no clear improvement is ob- 
tained but no information is lost, which again 
shows that there is no significant bias in the 
use of time as a third dimension. 

Figure 1 shows the 3D view of a simulated 
close- by event. The different volume shapes 
with dimensions coming from pixel and time 
bin size are clearly seen, the color codes show 
the time bin of detection associated to each 
volume. 

Energy and Profiles Reconstruction 

The 3D volumes in space correspond to the re- 
gions from which the observed photons were 
emitted. It is then possible to compare their 
charges to that given by a hypothesis - a lon- 
gitudinal and a lateral profile - taking into ac- 
count the light attenuation from the emission 



point (P) to the detector and the effective col- 
lection area seen by the telescope, as shown in 
equation 2. 

Expp = GH{-Kp, Xmax,Xo, A, N^a^) ■ (2) 
GO(Xp,Rp,XTOaa;) ■ RS{rp) ■ 

^ (47rrp2) 

where Yf is the fluorescence yield, GH is 
the longitudinal profile depending on the slant 
depth, Xp, and GO is the lateral profile, 
that depends also on Rp the distance to the 
axis (the profiles are given, respectively, by a 
Gaisser-Hillas and a Gora function [2], for ex- 
ample), RS is the Rayleigh Scattering factor, 
depending on the distance between the emis- 
sion and detection points rp (no Mie scattering 
is considered for now and multiple scattering 
contributions are neglected), and A represents 
the effective diaphragm area and may include 
constant calibration factors. 

In addition to the fluorescence light, also direct 
and scattered Cherenkov are included. Even 
for showers not directed to the eye, scattered 
Cherenkov can represent a non-negligible frac- 
tion of the detected light. The parameterisa- 
tion of the angular distribution function given 
in ref [5] is used for the direct contribution, 
while for the scattered one the convolution of 
the Cherenkov light production with Rayleigh 
scattering is considered. 

The light observed in the telescope is spread by 
the non-perfect optics in the mirror, producing 
a spot, which depends on the incidence angle 
and can be parameterised from simulation, in- 
cluding also camera shadow effects. This infor- 
mation can be included to find the exact posi- 
tions of the photons in the camera. Since the 
calibration is done by illuminating equally all 
the camera points, a factor can be found 
to c;orrect up the sensitive areas, and correct 
down the Mercedes regions. 

The Monte Carlo integration starts by finding 
a cylinder with the reconstructed shower axis, 
and enclosing the several reconstructed vol- 
umes. The expected number of photons from 
function 2 is evaluated for each of A''i x N^oi 
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randomly generated points, to cover all the 
Njjoi volumes found for each telescope. After- 
wards, N2 points are spread according to the 
spot parametrisation and, in each, the function 
is corrected by /,„ . The sum of the function for 
the final points in each volume, normalised by 
volume of the cylinder Vcyi , gives the expected 
number of photons per pixel and time bin, to 
be compared to the observed one: 
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Figure 2: Number of observed photons from a 
simulated event of 10^^-^ eV, as a function of 
time. The data points are compared to the to- 
tal expectations. The Cherenkov contribution 
is shown separately in the dashed line. 

Figure 2 shows the expected and observed 
number of detected photons as a function of 
the observed time, for a given simulated event. 
The simulated values are used for the Gaisser- 
Hillas and Gora function. Knowing the error 
on the observed value, one can create a func- 
tion to minimise for the model parameters ap- 
pearing in both functions. The lateral profile 
function could be directly tested for close-by 
events and that might help in the estimation 
of the total energy and Xmax , even when X^ax 
is out of the field of view. 



